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ABSTRACT 

We present an expanded sample of low-mass black holes (BHs) found in galactic nuclei. Using standard virial 
mass techniques to estimate BH masses, we select from the Fourth Data Release of the Sloan Digital Sky Survey 
all broad-line active galaxies with masses < 2 x 10* Mq. BHs in this mass regime provide unique tests of the 
relationship between BHs and galaxies, since their late-type galaxy hosts do not necessarily contain classical 
bulges. Furthermore, they provide observational analogs of primordial seed BHs and are expected, when merging, 
to provide strong gravitational signals for future detectors such as LISA. From our preliminary sample of 19, 
we have increased the total sample by an order of magnitude to 174, as well as an additional 55 (less secure) 
candidates. The sample has a median BH mass of (Mbh) = 1.3 x 10* Mq, and in general the objects are radiating 
at high fractions of their Eddington limits. We investigate the broad spectral properties of the sample; 55 are 
detected by ROSAT, with soft X-ray luminosities in the range lO'*" to 7 x 10"*^ ergs s"' . Much like the preliminary 
sample, these objects are predominantly radio-quiet (R = /6cm//4400 < 10), but 1 1 objects are detected at 20 cm, 
with radio powers (10^' - 10^^ W Hz"') that may arise from either star formation or nuclear activity; only 1% 
of the sample is radio-loud. We further confirm that, with {Mg) = -19.3 and {g- r) = 0.7 mag, the host galaxies 
are low-mass, late-type systems. At least 40% show disk-like morphologies, and the combination of host galaxy 
colors and higher-order Balmer absorption lines indicate intermediate-age stellar populations in a subset of the 
sample. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert 



1. INTRODUCTION 

There are many strong observational connections between 
galaxy bulges and central supermassive black holes (BHs). BH 
masses correlate surprisingly tightly with bulge properties, in- 
cluding luminosity (e.g., Marconi & Hunt 2003) and stellar ve- 
locity dispersion (Mbh - ; Ferrarese & Merritt 2000; Geb- 
hardt et al. 2000a; Tremaine et al. 2002). Also active galactic 
nuclei (AGNs) in the local Universe are predominantly found 
in massive, bulge-dominated galaxies (Ho et al. 1997b; Kauff- 
mann et al. 2003). Perhaps the formation of BHs and bulges 
are related, in which case we might not expect bulgeless galax- 
ies to obey scaling relations between BHs and galaxies, or even 
necessarily host a central BH. Indeed, dynamical study of the 
bulgeless spiral galaxy M33 places strong limits on the pres- 
ence of a dark central massive object (Gebhardt et al. 2001), 
while any BH in the nucleus of the dwarf spheroidal galaxy 
NGC 205 has been shown to lie below the low-mass extrap- 
olation of the Mbh - cr* relation (Valluri et al. 2005). On the 
other hand, the M31 globular cluster Gl shows dynamical and 
radiative evidence for a central BH (Gebhardt et al. 2002, 2005; 
Pooley & Rappaport 2006; Ulvestad et al. 2007). However, the 
mixed stellar population and high degree of rotational support 
in this massive cluster suggest that Gl is actually the nucleus of 
a tidally stripped dwarf galaxy (e.g., Meylan et al. 2001). Thus, 
dynamical studies present an ambiguous verdict on the pres- 
ence of nuclear BHs in dwarf stellar systems. Unfortunately, 
it is not currently feasible to spatially resolve the gravitational 



sphere of influence of a ^ 10^ Mq BH outside the Local Group, 
in order to search for BHs in more dwarf stellar systems using 
dynamical methods. 

Although they are difficult to find, the occupation fraction 
of nuclear BHs in dwarf systems and the space density of low- 
mass BHs are of considerable interest. Apart from furnishing 
additional insight into the possible origin of the Mbh - cr* re- 
lation, low-mass BHs provide low-redshift counterparts to the 
primordial seed BHs; the low-mass cut-off in the BH mass 
function today provides a constraint on the mass function of 
seed BHs. The merging of BHs in this mass range is expected 
to provide a strong signal for the gravitational wave experiment 
LISA (e.g., Hughes 2002). Furthermore, gravitational radiation 
recoil is expected to impart velocities to BH merger remnants 
that exceed the escape velocities of dwarf galaxies (e.g., Fa- 
vata et al. 2004; Merritt et al. 2004). This effect alone might 
decrease the occupation fraction of BHs in dwarf galaxies. 

In the absence of detectable dynamical signatures, we are 
forced to rely on less direct evidence for the presence of nuclear 
BHs, namely AGN activity. In fact, there are two well-studied 
AGNs in dwarf galaxies: NGC 4395 is a bulgeless spiral galaxy 
(Filippenko & Sargent 1989), while POX 52 (Kunth et al. 1987) 
is a dwarf spheroidal galaxy (Barth et al. 2004fl BH mass 
determinations are far less certain in the absence of dynami- 
cal constraints, but a variety of techniques yield BH masses for 
each object of - 10^ Mq (Filippenko & Ho 2003; Barth et al. 
2004; Peterson et al. 2005). Greene & Ho (2004; GH hereafter) 
performed a systematic search for AGNs with 



'Hubble Fellow and Princeton-Carnegie Fellow. 

^Such objects are commonly referred to as dwarf elliptical galaxies in the literature. However, because their structure is quite different from that of classical 
elliptical galaxies, we prefer to refer to them as dwarf spheroidal systems. 
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Table 1. The SDSS Sample 



ID 


SDSS Name 


Flag 


z 


9 


g-r 


A, 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1 


SDSS JOOOlll. 15-100155.5 




0.0489 


18.6 


0.7 


0.21 


2 


SDSS J000308.47+154842.3 




0.118 


18.2 


0.5 


0.16 


3 


SDSS J000605.59-092007.0 


c 


0.0699 


17.9 


0.7 


0.34 



Note. — Col. (1): Identification number assigned in this paper. Col. 
(2); Official SDSS name. Col. (3); Objects that are beneath our detection 
threshold are indicated with c. Previous identifier, if relevant, is indicated. 
Col. (4): Redshift measured by the SDSS pipeline. Col. (5): Petrosian 
g magnitude. Objects with questionable Petrosian magnitudes (fiagged as 
many-petro or no-petro by the SDSS pipeline) are noted with *, and the model 
magnitudes are tabulated. Col. (6): Petrosian g — r color. Col. (7): Galactic 
extinction in the g baJid. Table 1 is available in its entirety via the link to 
the machine-readable table above. The following is included only as a guide 
to content and presentation. 

BH masses < 10^ Mq, using the First Data Release of the Sloan 
Digital Sky Survey (SDSS; York et al. 2000) and recovered 19 
objects in that preliminary search. Remarkably, follow-up spec- 
troscopy using ESI on Keck suggests that these sources, as well 
as NGC 4395, POX 52, and Gl, are consistent with the low- 
mass extrapolation of the Mbh - relation (Barth et al. 2005). 
At the same time, many of the host galaxies are late-type galax- 
ies without classical bulges (J. E. Greene, et al. in preparation). 
Clearly, more objects are needed to both confirm this prelim- 
inary result, and also to constrain the structural parameters of 
the host galaxies. To this end, we have repeated our analysis 
using the Fourth Data Release of the SDSS (DR4; Adelman- 
McCarthy et al. 2006), and present the updated sample here. 

Our goal is to study the properties of the lowest-mass BHs 
that are identifiable with the SDSS. The upper mass limit on 
this sample is somewhat arbitrary. Within the inactive sam- 
ple of galaxies with dynamical BH masses, the BH in M32, 
Mbh=('2 -5^.2 0) X 10'' Mq, has the lowest mass (Tremaine et 
al. 2002; excluding the Milky Way, whose modem mass is 
3.5 X 10* Mq; e.g., Ghez et al. 2005; Eisenhauer et al. 2005). 
We have adopted 2x10*" Mq as the upper limit for our search. 
Because dynamical methods cannot extend into this regime, we 
previously had constraints on neither the existence of lower- 
mass BHs nor on whether they obey similar scaling relations 
with galaxy bulge properties as their high-mass cousins. Fur- 
thermore, while there are good reasons to believe that not every 
dwarf stellar system contains a BH, we do not know at what 
mass the occupation fraction departs from unity. 

Throughout we assume the following cosmological param- 
eters to calculate distances: //(> = 100 h = 71 km s~^ Mpc~', 
= 0.27, and = 0.75 (Spergel et al. 2003). 

2. SAMPLE SELECTION AND ANALYSIS 

Our low-mass BHs are selected from the sample of broad- 
line AGNs with z < 0.352 described in detail by Greene & 
Ho (2007b) and briefly reviewed here for completeness. We 
begin with DR4 of the SDSS and search for all AGNs with 
broad Ha, where "broad" in this context means a significant 
extra component relative to the narrow-line profile (based on 
the [S II] AA6716, 6731 doublet; e.g.. Ho et al. 1997c). Con- 
tinuum subtraction, which is crucial to uncover low-contrast 
broad lines, is performed with the principle component analysis 
method developed by Hao et al. (2005). Our subsequent selec- 



tion technique comprises a two-step procedure: we first select 
objects with high root-mean- square (rms) deviations above the 
continuum in the region potentially containing broad Ha, and 
then perform more detailed profile fitting to isolate those ob- 
jects with broad Ha profiles. Unfortunately, this selection pro- 
cess results in many sources of such low Ha luminosity that (1) 
their nature is ambiguous and (2) our ability to recover a reU- 
able BH mass is severely compromised. Based on simulations, 
we thus impose a combined rms-weighted flux and equivalent 
width (EW) cut designed to minimize spurious detections. The 
resulting sample comprises 8435 objects. 

Since we cannot use stellar dynamical methods to measure 
BH masses, we instead use the photoionized broad-line region 
(BLR) gas as a dynamical tracer of the BH mass. While the 
BLR velocity dispersion is derived from the line width, the 
BLR radius is inferred from the AGN luminosity (in this case 
Ha; Greene & Ho 2005b). The so-called radius-luminosity re- 
lation is derived from reverberation mapping of ~ 30 AGNs, 
for which radii are measured based on the delay between varia- 
tions in the AGN photoionizing continuum and BLR line emis- 
sion (e.g., Kaspi et al. 2005; Greene & Ho 2005b; Bentz et al. 
2006). The "virial" BH mass is simply MBH=fRv^/G, where 
/ is a scaling factor that accounts for the unknown geometry of 
the BLR, assumed here to be spherical (/=0.75; Netzer 1990). 
Although BH masses derived from reverberation mapping are 
susceptible to large systematic errors (due to uncertainties in 
the BLR geometry and kinematics; e.g., Krolik 2001), remark- 
ably they have been shown to agree with cr,, in those cases 
for which a* can be measured. The measured scatter is ~ 0.3 
dex (Gebhardt et al. 2000b; Ferrarese et al. 2001; Onken et al. 
2004; Nelson et al. 2004). Virial masses are less direct, since 
they rely on the radius-luminosity relation; in the largest such 
comparison to date, Greene & Ho (2006) find a scatter of 0.4 
dex of single-epoch virial masses about the Mbh - cr* relation. 
Our BH mass estimator is discussed in detail in the Appendix. 
Here we simply note that Greene & Ho (2007b) used an old 
version of the radius-luminosity relation slope, while here we 
have updated the value to that found by Bentz et al. (2006). In 
general, the derived masses increase by ~ 0.3 dex. 

Selecting those broad-line AGNs with masses < 2 x 10* Mq 
results in a sample of 174 objects, which are the main subject of 
this paper (Fig. 1). In addition, some (55) of the objects below 
the detection threshold are nevertheless interesting candidates. 
These objects have been selected by manual inspection of the 
spectra, and thus comprise a biased and incomplete sample. We 
worry, as explained above, lhal our BH masses are not stable 
for these Ha luminosities, but we deem these objects worthy of 
follow-up spectroscopy. In what follows we will refer to these 
objects as the c ("candidate") sample, and we will present all 
trends including and excluding these objects. Basic properties 
of the sample are suimnarized in Table 1. 

2.1. Fe II Fitting 

Our modeling of the H/3 region warrants some additional 
discussion. While the continua of our spectra, in general, 
are dominated by stellar light, in some cases there is an ad- 
ditional "pseudo-continuum" component contributed by broad 
Fe II multiplets. The Fe II emission extends over the entire opti- 
cal and UV spectrum, but is especially troublesome in the wave- 
length ranges 4400^800 A and 5150-5500 A. In principle, 
both the H/3 and [O m] AA4959, 5007 fits can be severely com- 
promised when the Fe II component is ignored (see Fig. 2). We 
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Fig. 1 . — Distributions of "virial" BH mass and Eddington ratio for tlie entire sample. Mbh is calculated from Ljja and FWHMna, using the formalism of Greene 
& Ho 2005b. The Eddington ratio is derived assuming an average bolometric correction of L|,oi = 2.34 X lO^CLHa/lO^^)"'** ergs s~' (see text). Filled histograms 
represent objects below the detection threshold defined in Greene & Ho 2007b (the c sample). 



follow Standard procedure and model the Fe II emission with 
an empirical template (e.g., Boroson & Green 1992; Greene & 
Ho 2005b). Since many of our sources also require substantial 
galaxy subtraction, the Fe II, H/3, and [O IIIJ hues are modeled 
simultaneously in the galaxy-continuum-subtracted spectrum, 
as described below. 

In addition to an Fe II template, our model includes narrow 
[O III], and narrow and broad H/3. In the standard way, all rel- 
ative wavelengths for central narrow components are fixed to 
laboratory values, and the narrow H/3 flux is fixed to be no more 
than 1/3.1 of the narrow component of Ha (Case B' recom- 
bination; Halpem & Steiner 1983; Gaskell & Ferland 1984). 
The [O III] lines are fit simultaneously with a core and a wing 
component, following Greene & Ho (2005a), with a relative 
strength of 2.96, while the narrow H/3 line is modeled with a 
single Gaussian due to its typical low signal-to-noise ratio. 

Usually, Fe II templates are derived from 1 Zw I, a strong 
Fe II source with relatively narrow broad lines. With a broad- 
line width of FWHM = 1240 km s'^Boroson & Green 1992), 
1 Zw I is significantly broader than many of the objects in our 
sample. For this reason, we are compelled to build our own Fe II 
template from the SDSS data. Using a high S/N spectrum of 
SDSS J155909.62+350147.4, whose FWHMh„ is 860 km s"', 
we model H/3+[0 III] using the 4750-5100 A region with no 
Fe II template included (Fig. 2a), and subtract this model. We 
further mask the regions 4840-4870 A and 5010-5018 A (note 
these wavelengths are in vacuum), since these regions are dom- 
inated by noise in the subtraction. While the resulting template 
is similar to 1 Zw I (Fig. 2b), it is clear that we do not re- 
cover the shape of the 5000 A feature perfectly, and that there 
are slight differences in continuum slope between the two tem- 
plates. To explore the importance of the former problem, we 
have repeated our fitting procedure using templates made with 
various masking regions, and find no difference in the derived 



fit parameters. The continuum shape, on the other hand, re- 
mains a systematic uncertainty inherent to our method; there 
is some degeneracy between the shape of the continuum and 
the amplitude of Fe II contamination, leading to uncertainties 
in the overall amplitude of Fe II. Nevertheless, our measured 
H/3h-[0 III] fits are robust to this uncertainty. When we fix the 
Fe n template to values ranging from one-half to twice the best- 
fit value, we find < 1% changes in all other measured param- 
eters. Presumably the [O III] is sufficiently resolved from the 
Fe II component at 5000 A so that small errors in the latter do 
not strongly impact our model of the former. 

2.2. Comparison with Greene & Ho (2004) 

Let us examine how the original GH sample compares with 
this new larger sample. Objects from the original paper will 
be referenced using their identification number in that paper 
(e.g., GH02). There are no significant differences between the 
two works, but there are some issues worth noting. In general, 
we have improved our selection procedure (see Greene & Ho 
2007b), so that our sample from DRl is close to two and a half 
times larger than the original sample. We are now using the Ha 
rather than the L5100 luminosity to derive the BLR radius, and 
we are using the updated radius-luminosity relation slope of 
Bentz et al. (2006). Only one (GH19) of the GH objects is not 
included in the final DR4 broad-line AGN sample of Greene & 
Ho (2007b), although GH07, GHll, GH15, GH16, and GH18 
now have masses above our mass boundary of 2 x 10*" Mq. Fi- 
nally, the original GH sample excluded galaxies with significant 
galaxy contamination, due to the fear that these objects would 
comprise more massive BHs with low-contrast broad lines that 
were rendered undetectable by the large galaxy luminosity. In 
this work we have not explicitly included such a cut, but our 
EW threshold is operationally similar. While some of the ob- 
jects in the c sample may consist of high-mass interlopers, 
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Fig. 2.— (a) The H/3+[0 iii] fit to SDSS J155909.62+350147 used to create 
the Fe ii template shown in (h). Shown are the continuum-subtracted (top), 
best-fit {red thin) and residual (bottom) spectra. Units are 10"'^ ergs s"' cm"^ 
A"', (b) A compai'ison of the standard Boroson & Green 1992 1 Zw i Fe ii 
template (blue dotted) with that used in this work (solid). Spectra have been 
normalized to unity. The two templates are very similar, although our template 
is clearly narrower, and, as discussed in the text, the overall slope is somewhat 
different, (c) A sample fit to the H/3+[0 iii] region. Spectra in solid are the 
continuum-subtracted (top) and residual (bottom) spectra. The total model is 
oveiplotted (red thin), along with individual components (blue dotted). Units 
the same as (a) above. 

we find in general very good agreement in properties between 
the two samples. 

2.3. Comparison with Dong et al. (2007a, b) 

Dong et al. (2007a) have presented a low-mass BH with a 
mass of 7 X lO'* Mq. This source was selected from the Fifth 
Data Release of the SDSS, and so cannot be cross-checked with 
our sample. However, the object has been serendipitously ob- 
served with HST and XMM-Newton, as well as being detected 
with the Roentgen Satellite (ROSAT). There is tentative evi- 
dence for variability in the X-ray source (Lx = 7 x 10"^*^ ergs s~') 



and the host galaxy is a disk galaxy with Mr = -17.8 mag and a 
nuclear bar. The same group has identified a sample of 245 low- 
mass BHs from DR4, with a mass range of 5 x 10"*- 10'' Mq, 
and Eddington ratios ranging from 0.02 to 8, comparable to our 
sample (Dong et al. 2007b). Samples selected from two inde- 
pendent methods will provide important tests on the limitations 
of our method. Note, however, that the published masses are 
based on the radius-luminosity relation slope of Kaspi et al. 
(2005), and would increase by ~ 0.3 dex according to the for- 
malism adopted in this paper. 

3. SAMPLE PROPERTIES 

The final sample has a median redshift of 0.086 (0.099 if 
the c sample is excluded). The original sample has a slightly 
lower median z = 0.08 due in part to subsequent deep obser- 
vations in the Southern Strip (targeting quasars to a limiting 
magnitude of / = 19.9 mag rather than / = 19.1 mag). In Figure 
1 we show the distribution of "virial" BH mass and Edding- 
ton ratio for the sample, where L^dd = 1.26 x 10^** (Mbh/Mq) 
ergs s"'. Assuming that Lboi = 9.8 Lsioo (McLure & Dun- 
lop 2004), in terms of Lea our bolometric correction is Lboi = 
2.34 X lQ*\LHa/l0^y-^^ ergs s"' (Greene & Ho 2005b). The 
median BH mass is (Mbh) = 1-3 x 10^ Mq. With a median 
(Lboi/iEdd) = 0.4, we are clearly dominated by sources radiat- 
ing at substantial fractions of their Eddington limits (Fig. 1). 
This is not surprising; as discussed by Greene & Ho (2007b), 
we are only sensitive to the most luminous BHs in this mass 
regime. 

3.1. Optical Spectral Properties 

The statistical power of our new enlarged sample provides 
new constraints on the ensemble physical properties of both the 
radiating BHs and their host galaxies. We begin by placing 
the sources on the two-dimensional diagnostic diagrams that in 
combination discriminate between a stellar ionizing spectrum 
or one considerably harder (e.g., Baldwin et al. 1981; Veilleux 
& Osterbrock 1987). These line ratios are typically used to di- 
vide emission-line galaxies into H II galaxies and narrow-line 
AGNs (e.g.. Ho et al. 1997a; Kauffmann et al. 2003; Hao et 
al. 2005; Table 2). In contrast, broad-line-selected objects are 
unbiased with regard to position on the diagnostic diagrams, 
while at the same time we know that some fraction of the ion- 
izing continuum is contributed by an AGN. 

We focus first on the [N II] A6583/Ha diagnostic diagram 
(Fig. 3a). We do not see many objects in the lower-right quad- 
rant. Objects there are low-ionization nuclear emission-line re- 
gion (LINER; Heckman 1980) sources and typically are highly 
sub-Eddington AGNs (Ho 2004). Compared to classical Seyfert 
galaxies, which occupy the upper-right region of the diagnostic 
diagram, we see that our sources span a broader range in both 
[N II]/Ha and [O III]/H/3. The decrease in [N ll]/Ha is most eas- 
ily explained as a decrease in the gas-phase metallicity of these 
AGNs (e.g., Kraemer et al. 1999; Groves et al. 2006). Since ni- 
trogen is a secondary element, [N II] is particularly sensitive to 
changes in metallicity for > O.IZq. In contrast, the [O III]/H/3 
ratio changes much more slowly, because as the metallicity 
decreases, the temperature of the NLR increases correspond- 
ingly, which tends to boost the [O III] strength (e.g.. Groves 
et al. 2004a, 2006; see metallicity tracks on Fig. 4). Rela- 
tively low metallicities for this sample are not unexpected, since 
local AGN hosts in general are known to be massive, bulge- 
dominated galaxies (e.g.. Ho et al. 1997b; Kauffmann et al. 
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5 



ID 


[O 11] 


Fe II 


(H0)„ 


(H/3)b 


[O III] 


[Or] 


Ha„ 


Hc«b 


[Nil] 


[S II] 


[S II] 


FWHMhc 


FWHM[oiii] 




A3727 


A4570 






A5007 


A6300 






A6583 


A6716 


A6731 


(km s-i) 


(km s-i) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


1 


0.46 


0.18 


0.27 


<0.41 


1.84 


0.067 


1.4 


1.3 


0.51 


0.31 


0.20 


1870 


230 


2 


1.2 


1.6 


0.57 


1.0 


0.971 


0.16 


1.9 


4.8 


0.79 


0.43 


0.33 


976 


162 


3 


0.53 


<0.39 


0.18 


0.62 


1.88 


0.10 


0.69 


0.84 


0.58 


0.24 


0.20 


1710 


215 



Note. Col. (1): Identification number assigned in tliis paper. Col. (2) (12): All fluxes are relative to that of [O III] A5007, which 
is listed in units of lO^^^ergs s^^ cm^^. Note that these are observed values; only Galactic extinction correction has been applied. The 
subscripts "n" and ''b" in Col. (4) (5) and (8) (9) refer to the narrow and broad components of the lino, respectively. Col. (13) -(14): 
Linewidths have been corrected for instrumental resolution using the values measured from arc spectra and tabulated by the SDSS. 
Table 2 is available in its entirety via the link to the machine-readable table above. The following is included only as a guide to content 
and presentation. 



2003) , while these low-mass BHs tend to live in low-mass and 
correspondingly low-metallicity systems (§3.4; Tremonti et al. 

2004) . 

More striking is the apparent corresponding decrease in 
[O III]/H/3, which as we have seen, cannot be ascribed to a de- 
crease in metallicity. In fact, |0 III]/H/3 appears to decrease in 
concert with [N II]/Ha, [S II]/Hq, and [O I] A6300/Ha (Fig. 
3b, c). At the same time, there is a clear trend of increasing 
[O IIJ A3727/[0 III] line ratios as the [O III|/H/3 ratio decreases 
(Fig. 4). Changes in a variety of NLR conditions, including 
metalhcity, electron density (rie), ionization parameter (U), ion- 
izing spectral shape, and reddening might be invoked to explain 
the observed line ratios, while particular properties of either the 
AGN (e.g., luminosity or Eddington ratio) or the host galaxy 
(e.g., star formation) may be responsible for the changing con- 
ditions. As we will argue, it is unlikely that a single parameter 
is responsible for all of these trends. Correlated errors might 
also be expected to move line ratios in particular ways across 
the diagnostic diagrams. If, for instance, we are systematically 
overestimating the strength of narrow Ha in systems with weak 
NLR emission, then we might expect correlated tracks as ob- 
served. However, [O III] is never blended with broad H/? in 
these sources, and thus [O III] errors are not correlated with 
the other lines, but nevertheless the [O III] line strength spans a 
wide range relative to H/3. 

One might imagine that the intrinsic Une ratios are all lo- 
cated in the Seyfert locus, but that contamination from star for- 
mation spreads out their positions in the diagnostic diagrams. 
We perform a simple thought experiment to demonstrate that 
there must be an intrinsic spread in line ratios to reproduce 
the observed trends. Since star-forming galaxies follow a well- 
defined sequence in the diagnostic diagram (bounded by the 
Kauffmann et al. 2003 relation shown Fig. 5), we can eas- 
ily investigate the type of bias that star formation would con- 
tribute. In Figure 5 we have chosen three fiducial AGN posi- 
tions in the diagnostic diagram. Position A represents a typ- 
ical high-ionization system, while B has low [O III]/H/3 and 
[N n]/Ha and C represents a LINER. To each AGN we add 
varying amounts of star formation with fine ratios drawn from 
the Kauffmann line, and the H/3 strength in star formation vary- 
ing from half to twice the AGN H/3 luminosity. The tracks 
(shown as dash-dot, solid, and dashed lines for positions A, 
B, and C, respectively) demonstrate that while star formation 
may increase the spread in positions on the diagnostic diagram, 
it alone cannot move objects from the high-ionization position 
A to the low-[0 IIIJ/H/3 positions we observe. There must be a 



true spread in the intrinsic AGN line ratios. Furthermore, posi- 
tion C is disfavored, since only with a large fraction of the light 
contributed by low-metalUcity star formation might we move 
the objects from position C toward position B. If the objects in 
this sample follow the mass-metallicity relation of Tremonti et 
al. (2004), with typical luminosities of Mg « -19.3 (Table 3, 
§3.4), they would not, on average, have sub-solar metallicity. 

Given that the dispersion in line ratios must be intrinsic to 
the NLR, we investigate various means to explain the observed 
line-ratio distributions. Ho et al. (1993) find that narrow-line 
AGNs display a similar trend in the occupation of diagnostic 
diagrams, which they ascribe to a trend of harder ionizing spec- 
tral shape at lower luminosity. Exactly this trend has been seen 
in the X-ray to optical slopes of broad-line AGNs over a wide 
range in luminosity (e.g., Strateva et al. 2005), including for a 
subset of the GH sample (Greene & Ho 2007a). A correlation 
between low-luminosities and harder spectral shape would lead 
to a correlation between AGN luminosity and line ratios, but 
we see no trends between AGN luminosity (as traced by Lua) 
and either [O III]/H/3 or [O I]/Ha. 

We have already discussed the impact of varying metalhc- 
ity above, and since [N II]/Ha is most sensitive to its varia- 
tion, metallicity is unlikely to explain the correlated decrease 
in [S Ill/Ha, [O I]/Hq! and [N II]/Hq:. Furthermore, metallic- 
ity variations alone cannot explain the large observed range 
in [O n]/[0 III] ratios (Fig. 4). Changes in either ioniza- 
tion parameter or density might explain the [O n]/[0 III] ra- 
tios; increases in either leads to a relative increase in [O III]. 
Indeed, electron densities inferred from the ratio of the [S II] 
Unes show a significant correlation with [O III]/H/3. The Spear- 
man rank correlation coefficient is p = 0.296, with a probabil- 
ity P < 10"^ that no correlation is present. However, there is 
only a mildly significant trend between electron density and 
[N II]/Ha (p = 0.197, P = 0.003) and no trend with [S n]/Ha 
ip = 0.124, P = 0.06; as one might expect given their lower 
critical densities). On the other hand, changes in ionization pa- 
rameter move, for example, [N II]/Ha vs. |0 III]/H/3 on tracks 
perpendicular to those observed (e.g., Groves et al. 2004a, b). 
A single parameter seems unable to reproduce all of the ob- 
served trends simultaneously. Presumably this is due, in part, 
to stratification in the density and ionization parameter of the 
NLR with radius, which leads to different typical conditions 
in the emission regions of different lines (e.g., FiUppenko & 
Halpern 1984). 

In summary, a single parameter cannot be invoked to explain 
the observed positions in the diagnostic diagrams. 
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Fig. 3.— Diagnostic diagrams plotting log [O iii] A5007/H/3 versus (a) log [N ii] A6583/Ha, (b) log [S ii] AA6716, 6731/Ha, and (c) log [O i] A630O/Ha. The 
c sample is shown in open symbols, and there are red boxes around ROSAT detections and blue triangles around FIRST detections. The line ratios have not been 
corrected for reddening, but this should not matter because of the small wavelength separation of the lines. The dotted lines mark the boundaries of the three main 
classes of emission-line nuclei, according to the convention of Ho et al. 1997a; "H" = H ii nuclei, "S" = Seyferts, and "L" = LINERs. 



Table 3. Luminosity and Mass Measurements 



ID (total) 


A/g(AGN) 


A/9{host) 


log Lhc 


log il/BH 


log .^bol/i-Edd 


log Mold 


(1) (2) 


(3) 


(4) 


(5) 


(8) 


(7) 


(8) 


1 -18.02 


-15.91 


-17.82 


40.11 


6.2 


-1.6 


5.8 


2 -20.16 


-18.24 


-19.91 


41.22 


6.1 


-0.5 


5.9 


3 -19.44 


-16.31 


-19.37 


40.26 


6.2 


-1.4 




Note. — Col. (1) 


Identificatic 


n number 


assigned in 


this paper. 


Col. (2): Total g 


-band absolute 



magnitude. Col. (3): AGN 5-band absolute magnitude, estimated Irom Lh^ given in Col. (5) 
and a conversion from L^a to Mg assuming fx oc A^^'^^. Col. (4): Host galaxy p-band absolute 
magnitude, obtained by subtracting the AGN luminosity from the total luminosity. Cases where 
the AGN accounts for the total galaxy luminosity are indicated with an ellipsis. Col. (5): AGN 
luminosity in broad Ka (ergs s-^). Col. (6): Virial mass estimate of the BH (Mq). Col. (7): Ratio 
of the bolometric luminosity (see text) to the Eddington luminosity. Col. (8): Previously reported 
virial mass estimate of the BH (Mq) from Greene & Ho 2007b. There are no entries for the c sample. 
Table 3 is available in its entirety via the link to the machine-readable table above. The following is 
included only as a guide to content and presentation. 

Although changing physical conditions (such as electron den- 
sity) within the NLR may be important, the objects lying close 
to the star-forming locus may simply have an extremely weak 
NLR, such that the majority of the narrow emission is com- 
ing from H II regions. A similar phenomenon, the "vanishing" 
NLR, has been seen in a substantial fraction of higher redshift 
AGNs at high Eddington ratios (e.g., Netzer et al. 2004, 2006). 
If observed correlations between NLR size and AGN luminos- 
ity (e.g., Bennert et al. 2002; Schmitt et al. 2003) are extrap- 
olated to the highest luminosities, then one would predict un- 
bound nebulosities; Netzer et al. (2004) suppose that the NLRs 
in these systems grew so large that they were no longer bound 
to their host galaxies. The luminosities in our systems are sig- 
nificantly lower, and the expected NLR sizes are < 100 pc if 
they obey the NLR size-luminosity relation, but we cannot rule 
out that the AGN actually expelled some large fraction of the 
surrounding ISM, leading to overall weakness in the line inten- 
sities. 

Finally, there is an apparent connection between [O III]/H/3 
and Eddington ratio (e.g., Boroson 2002). However, for 
these objects there is no trend between [O III]/H/3(narrow) and 
Lboi/LEdd ip = 0.043, P = 0.5) or [O III] equivalent width and 
i'boi/iEdd ip = 0.093, P = 0.2). These trends do not change 
when the c sample is excluded. On the other hand, although 
no physically motivated explanation exists, low [O III]/H/? tra- 
ditionally was a defining characteristic of the class of objects 
known as narrow-line Seyfert 1 galaxies (NLSls; Osterbrock 
& Pogge 1985), 



Table 4. ROSAT Detections 



ID 


C 


log Nn 


log fx 


log Lx 


(1) 


(2) 


(3) 


(4) 


(5) 


18 


0.023± 0.011 


20.55 


-12.66±0.16 


42.44±0.16 


24 


0.053± 0.019 


20.86 


-12.17±0.05 


42.49±0.05 


28 


0.024± 0.0090 


20.79 


-12.54±0.07 


42.19±0.07 



Note. — Col. (1): Identification number assigned in this pa- 
per. Col. (2): R.OSAT count rate (counts sec^^). Col. (3): 
Galactic column A^h (cni^^), calculated following Dickey & Lock- 
man 1990. Col. (4): X-ray flux in the 0.5-2 keV band (ergs s^^ 
cm^^) assuming a power-law spectrum with Ts — 1 and A^h 
from col. (3). Col. (5): X-ray luminosity in the 0.5—2 keV band 
(ergs s^^). Table 4 is available in its entirety via the link to the 
machine-readable table above. The following is included only as 
a guide to content and presentation. 

which are broad-line AGNs with relatively narrow broad-line 
widths. We turn now to a general comparison of our objects 
with NLSls. 

3.2. Comparison with NLSls 

While NLSls are an observationally defined class, a pic- 
ture has emerged, based on their strong Fe II/H/3 and weak 
[O III]/H/3 line ratios, strong soft X-ray excess (e.g.. Boiler et al. 
1996) and X-ray vai-iabihty (Leighly 1999), that NLSls are in 
general low-mass BHs (thus accounting for the relatively nar- 
row broad lines) in a high accretion state (e.g.. Pounds et al. 
1995). Since the current sample is selected based on BH mass 
(and indirectly Lboi/i^Edd), it represents a very uniform, opti- 
cally selected NLS 1 sample. In this section, we highlight com- 
parisons between the optical and broad spectral properties of 
this and other NLS 1 samples. 

One of the supporting pieces of evidence that NLSls are ra- 
diating close to their Eddington limit comes from their observed 
high Fe II/H/3 ratios (e.g., Veron-Cetty et al. 2001; throughout 
this section we refer to total [broad-i-narrow] H/3 luminosity). 
Using principle component analysis, Boroson & Green (1992) 
found that NLSls lie at one extreme of their first component 
(fondly known as Eigenvector 1), and later work has shown 
that Eigenvector 1 properties, including high Fe II/H/3 and low 
[O III]/H/3 ratios. 
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Fig. 4. — Tracks in [O ii]/[0 iii] for different metallicities, fol- 
lowing Groves et al. 2006. At a given metallicity, tlie tracks span a 
range in ionization parameter tJo = 5«/(ciiH), wliere S« is the number 
of ionizing photons, and iin = Pgis/iklO*). The model values are log 
Uo=(-0.40,-1.03, -1.43, -1.88, -2.24, -2.57, -2.98) and the tracks are labeled 
at the low-Uo end, as indicated with dashed lines. It is clear that metallicity 
alone cannot account for the spread in [O ii]/[0 iii] ratios we observe. Open 
symbols indicate the c sample. 

as well as soft X-ray excess and radio weakness, tend to be cor- 
related with Eddington ratio (e.g., Boroson 2002). On the other 
hand, GH found that their objects span a larger range in both 
Fe II and [O III] strengths relative to H/3 than classic NLSls. 

The current sample has somewhat lower Fe II/H/3 strengths 
than GH found, perhaps because there are objects with more 
significant galaxy contamination in this sample. Using the 
Kaplan-Meier product-limit estimator, which accounts for up- 
per limits (Feigelson & Nelson 1985), we find (Fe II/H/?) = 
0.57 ±0.03 (0.61 ±0.03) with (without) the c sample, as com- 
pared to (Fe II/H/3) = 0.67 ±0.04 for the 56 NLSls presented 
by Veron-Cetty et al. (2001), and (Fe II/H/3) = 1.32 ±0.16 for 
the original GH sample. The [O III] distributions are consistent 
with our previous findings: ([O III]/H/3) = 2.6 ±0.3 (1 .9 ±0.2), 
while we found ([O III]/H/3) = 2.24 ± 0.72 for the original sam- 
ple. Although it has been customary to cite [O III]/H/3(total), 
the ratio of [O III]/H/3(narrow) is more straightforward to inter- 
pret, and we report it as well: ([O III]/H/3(narrow)) = 6. 1 ± 0.5 
(6.5 ± 0.6). In all other discussion, we report trends for 
[O III]/H/3(narrow). 

More generally, trends in spectral properties with Eddington 
ratio may be manifest in composite spectra. We investigate four 
bins of increasing Eddington ratio (similar to increasing bins in 
luminosity over such a restricted range in BH mass; Fig. 6), 
and we have included the c objects, which comprise the lowest 
Lboi/iEdd bin. To compute the composite, each spectrum is nor- 
malized to the median value in the spectral regions 4100-4800 
A and 5100-5800 A, resampled onto the same wavelength grid 
and then median combined. Median stacking preserves rela- 
tive line ratios, but not the continuum shape (e.g., Francis et al. 
1991). Although we do not explicitly weight the spectra, they 
are each normalized to their respective continuum level, which 
gives extra weight to faint sources. Errors at each pixel are cal- 
culated as the 68% interquartile range of the included pixel 
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Fig. 5. — A diagnostic diagram highlighting how different effects move 
points in the parameter space. For pure AGNs with line ratios originating 
at positions A, B, and C, we have mixed star-forming galaxies with line ra- 
tios along the Kauffmann et al. 2003 line (blue dotted), assuming that the H/3 
strength from photoionization by starlight ranges from half (upper track; A 
in long dash-dot, B in solid, C in long dash) to twice (lower track) that by 
the AGN. The thin arrows indicate the trajectory of increasing contributions 
from star formation. Also shown are photoionization models from Groves et 
al. 2006. Each track covers a single metallicity and the same range of ioniza- 
tion parameters as Figure 4 above. While mixing stellar and AGN-excited line 
ratios may broaden the locus of points in the diagnostic diagrams, it is clear 
that the intrinsic line ratios do not all originate at position A. 

values, following, e.g.. Fine et al. (2006). Note that this pro- 
cedure is considerably simpler than that of, e.g., Francis et al. 
(1991), since the redshifts measured by SDSS are highly re- 
liable and the spectral coverage and depth are so uniform for 
the majority of the sample. The clearest trend as we move 
to higher Eddington ratio is the decreasing importance of the 
galaxy continuum. At the same time, the Fe II features are 
much more apparent in the higher Eddington ratio bins. This 
appears to be in keeping with our expectations from Eigenvec- 
tor 1. On the other hand, using the original, unbinned data, 
we do not find significant evidence for a correlation between 
Fe II/H/3 and Lboi/-LEdd ip = -0.07, P = 0.3). There is a de- 
crease in the [O III]/H/3(total) ratio with Eddington ratio in the 
composite spectra, which ranges from 3.6 ± 1.6 in the lowest 
Lboi/i'Edd bin to 0.65 ±0.03 (or 2. 4 ±0.3 when considering only 
the narrow H/3). Errors are dominated by the variance in line 
strengths for individual objects comprising the composite. 

3.3. Radio and X-ray Spectral Properties 

There is additional information about the accretion proper- 
ties of the sample in their broader spectral energy distributions. 
NLSls, in general, are known for their strong soft X-ray ex- 
cess, but Greene & Ho (2007a; see also Williams et al. 2004) 
have found that objects selected purely on the basis of mass are 
not as extreme in this regard as classic NLSls. In this case, us- 
ing 5 ks observations of the 10 nearest GH objects, they found 
that the objects had very typical 0.5-2 keV power-law slopes of 
Fj w 1.8 [where N{E) oc E~^']. Within the current sample, 55 
are detected by the ROSAT All-Sky Survey (Voges et al. 1999), 
of which only two are in the c sample, with fluxes ranging from 
5.9 X 10"'^ to 3.4 X 10"'^ ergs s"' cm"^. Fluxes are derived 
from the ROSAT count rates using WebPimms (Mukai 1993) 
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Fig. 6. — Composite spectra in bins of increasing Eddington ratio, as shown. Each spectrum is the median combination of all spectra in the indicated Lbol/^'Edd 
bin, with the total number of galaxies indicated. Each spectrum is normalized to emission-line-free bands prior to combination (see text for details). Note the 
increasing strength of the Fe ii features and decreasing galaxy continuum strength with increasing Eddington ratio. 



and assuming a slope of Fj = 2. In calculating uncertainties in 
the fluxes, Fj is allowed to vary from 1 to 3, which is the full 
range observed by Greene & Ho (2007a). The con^esponding 
range of 0.5-2 keV luminosities is 10"*" to 7 x lO"^-' ergs s"' (Ta- 
ble 4). Since inactive galaxies with optical luminosities similar 
to those of the galaxies in our sample (Lb ~ 10''- 10'" L©) are 
expected to have typical X-ray luminosities of lO"'^ ergs s"' 
(e.g., Fabbiano 1989), the AGN most likely dominates the X- 
ray emission from the majority of these sources. 

Customarily, the ratio of optical to X-ray luminosity is ex- 
pressed in terms of the slope of a supposed power-law between 
the optical and X-ray bands, ao^= -0.38381og(/25()o//2kev)- 
The flux density at 2500 A, /2500, is calculated using the mea- 
sured Lua, the relation between Lea and Lgioo from Greene 
& Ho (2005b), and an assumed spectral slope of /3 = -1.56 
(/a oc A"'^; Vanden Berk et al. 2001). For our sample, the 
values lie in the range -1.5 <aox< -0.69, with a median of 
(ttox) = -1 .04. We would expect an aox~ -1 -2 corresponding to 
the median /25{)() = 6 x 10^^ ergs s"' Hz"' (as was seen for the GH 
objects observed with Chandra; Greene & Ho 2007a). While 
the observed median is slightly higher than one would expect 
(e.g., Strateva et al. 2005; Steffen et al. 2006), it is within the 
observed scatter, and given that these are the ROSAT-detected 
members of our sample, it is not surprising to find that they 
are X-ray bright. We see no obvious trends between and 
Lboi/iEdd or Fe II/H/3, but we caution that the dynamic range of 
this sample is very limited, and that we do not have meaningful 
upper limits for the X-ray-undetected sources. 

Radio emission is a complementary probe of the accretion 
process, if difficult to interpret cleanly. In general, NLSls are 
radio-quiet (see review in Greene et al. 2006), where radio-quiet 
typically means R < 10 (R = /6cm//440o)- Weak radio emission 



is part of Eigenvector 1 (e.g., Boroson & Green 1992), and mir- 
rors the general tendency for BHs of all masses to be radio-quiet 
when in high accretion states (e.g., McClintock & Remillard 
2006). Also, Best et al. (2005) find a strong preference for ra- 
dio sources to live in massive galaxies. A low incidence of radio 
activity in this sample is thus unsurprising. Eleven of the ob- 
jects are detected by the VLA0 FIRST survey at 20 cm although 
all but 18 objects are within the FIRST footprint (Becker et al. 
1995. The radio sources at 20 cm are unresolved in all cases, 
with radio powers that range from 2 x 10^' to 9 x 10^^ W Hz"' 
(Table 5). Inactive star-forming galaxies can easily produce this 
level of radio emission (e.g., Condon 1989), which need not be 
AGN-dominated. The exception is GHIO, which we have im- 
aged with the VLA in A-configuration at 6 cm (Greene et al. 
2006). We find that the source remains completely unresolved, 
suggesting that the AGN provides the dominant source of emis- 
sion from this galaxy. 

Of the objects with radio emission, four are in the c sample 
and five others are radio-quiet, which translates into a radio- 
loud fraction of 1% (accounting for the 18 objects in the sam- 
ple that are outside of the FIRST footprint), to be compared 
with the range of 0%-6% for NLSls in the literature (Greene 
et al. 2006). In Figure 7 the detected sources are placed on the 
relation between radio power and [O III] luminosity observed 
for local AGNs. They are roughly consistent with the relation 
measured for Seyfert galaxies (Ho & Peng 2001), but they are 
systematically higher than the fiducial relation. In addition, we 
have retrieved FIRST cut-outs around each non-detection. We 
use these images to derive upper limits on the radio luminosity 
(Table 5), and we also stack the images in two bins divided at an 
[O III] luminosity of L[oiii]=3 x lO''" ergs s"'. Because FIRST 
is a uniform survey, no weighting is applied, but in 16 cases 
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Table 5. FIRST Detections 



ID 


'S'20cm 


log -P20cm 


log R 


(1) 


(2) 


(3) 


(4) 


47 


1.77±0.14 


21.41±0.04 


0.65 


69 


1.43±0.15 


22.87±0.05 


0.96 


87 


1 no 1 n eio 
L.[JO±U.2Z 




0.60 


101 


c\ At I n t A 




1.90 


106 


5.yD±0.17 


22.1y±0.01 


1.30 


140 


1.03±0.15 


21.94±0.07 


1.20 


146 


1.29±0.15 


22.31±0.05 


0.44 


158 


1.55±0.15 


21.87±0.04 


0.89 


163 


2.33±0.14 


21.37±0.03 


1.10 


174 


1.91±0.15 


22.95±0.04 


1.20 


203 


3.39±0.14 


21.87±0.02 


0.00 


1 


<0.33±0.16 


<21.27±0.23 


<0.71 


2 


<0.45±0.15 


<22.20±0.15 


<0.70 


3 


<0.20±0.11 


<21.37±0.28 


<0.69 



Note. — Col. (1): Identification number as- 
signed in this paper. Detections are followed by 
upper limits. Col. (2): Flux density (or upper 
limit) at 20 cm from FIRST (mjy). Col. (3): 
Corresponding radio power (W Hz^-^). Col. 
(4): R = /6cm//^^gp^, assuming a spectral in- 
dex of (ii- = 0.46 (radio; Ulvestad & Ho 2001) 
and Qo = 0.44 (optical: Vanden Berk et al. 
2001), where oc 1/"°". Table 6 is available in 
its entirety via the link to the machine-readable 
table above. The following is included only as 
a guide to content and presentation. 

there is a ^ 2 (t detection in the inner five square pixels (where a 
is measured in a box off the center of the cut-out). We exclude 
these from our stacked images, although the answer changes 
negUgibly when they are included. We do detect radio emis- 
sion in each stacked image, shown as crossed squares in Figure 
7, also consistent with the general relation. 



3.4. Host Galaxies 

Naively, our expectation is that lower-mass BHs will be 
found, in general, in low-mass galaxies. However, the tight ob- 
served relations between BHs and galaxies pertain specifically 
to galaxy bulges. We are now in a position to test whether BH- 
bulge relations break down in low-mass galaxies. The spectra 
provide an estimate for the AGN luminosity (we prefer to use 
the Ha luminosity because of its robustness in the presence of 
significant galaxy contamination; Greene & Ho 2005b), which 
we convert to magnitudes using SDSS filter functions and as- 
suming a spectral shape of /? = -1 .56 (Vanden Berk et al. 2001). 
If we assume a flatter slope of /? = -1, as suggested by GH, the 
median magnitude difference is only 0.06 mag. In addition to 
the power law, we also include the observed strong emission 
Unes (H/3, [O III], [N II], and Ha) in the calculation of the AGN 
luminosity (on average, line emission accounts for (g) « 0.08 
mag). Once we have removed our best estimate for the AGN lu- 
minosity from the SDSS Petrosian magnitudes, we calculate K- 
corrections using the routines described in Blanton & Roweis 
(2007). The distribution of host galaxy luminosities thus de- 
rived has a median host galaxy luminosity of (Mg) = -19.3 mag 
(with or without the c sample) and is shown in Figure 8. With 
this method of AGN removal, we tend to underestimate the 
AGN luminosity due to fiber losses, and thereby 
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Fig. 7. — Adopted from Greene et al. 2006, this figure shows the rela- 
tion between Z,[oiii] ™d radio power for a large sample of AGNs. Objects 
from this paper are shown as open stars, while stacked luminosities for unde- 
tected sources are shown as crossed squares and upper Umits from the origi- 
nal GH sample are shown as open triangles. For clarity, our data are plotted 
alone in the inset. The solid line represents the fit from Ho & Peng 2001 to 
radio-quiet Seyfert galaxies and Palomar-Green (PG; Schmidt & Green 1983) 
quasars: Peom = (0.46 ± 0.15) Z,[oni] + (2.68 ± 6.21). The Ho & Peng sample 
is overplotted with filled (radio-loud sources) and open (radio-quiet sources) 
symbols. PG quasars are shown as squares, and Seyferts are shown as circles. 
NGC4395, withL[oiii]=5.8x 10^' ergs s"' (adjusted to a distance of 4.2 Mpc; 
Thim et al. 2004) is highlighted as a semifiUed circle. 

(conservatively) overestimate the galaxy luminosity. As a lower 
bound on the galaxy luminosity, we use the SDSS point-spread 
function (PSF; unresolved) magnitude as an alternate estimate 
of the AGN luminosity, and then apply a .fiT-correction to the 
corrected colors in the same manner. This method, on aver- 
age, overestimates the AGN luminosity, and results in a median 
host galaxy luminosity of (Mg) = -19.0 mag (Fig. 8, dotted 
histogram). The true host luminosities lie somewhere between 
these two bounds. For reference, Blanton et al. (2003) find a 
characteristic luminosity of M* = -20. 1 mag at z = 0. 1 (for our 
assumed cosmology), and thus our targets are 0.8 - 1 . 1 mag 
fainter than M* . 

Morphological information is difficult to extract from the 
SDSS images alone, due to their limited depth and spatial res- 
olution, but it is clear that BHs selected to have low masses are 
found in low-mass galaxies. However, if we restrict our atten- 
tion to the 139 sources with z < 0.1, we find ~ 80 clear disk 
galaxies (see GH for example SDSS images for the original 
sample). We can also calculate galaxy colors: (g-r) = 0.7 mag 
for the spectroscopic estimate of the AGN luminosity (0.6 mag 
for the PSF method). According to Fukugita et al. (1995), the 
host galaxies have the colors of typical Sab galaxies (g-r = 0.66 
mag). HST images of the GH sample indicate that there are 
actually two different populations of host galaxies. One is 
composed of late-type spiral galaxies, typically with strong 
bars. The other — the majority — group is comprised of compact 
spheroids, not unlike POX 52, the prototypical dwarf spheroidal 
galaxy hosting an AGN, which has red colors (B-V = 0.8 or 
g - r « 0.67 mag). We do not yet know if the spheroidal galax- 
ies in the larger sample are red in general. 
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3 .5 . Stellar Populations and Ongoing Star Formation 

It is tempting to imagine that the constant mass ratio ob- 
served between BHs and bulges results from concurrent (or 
at least orchestrated) star formation and AGN activity. In- 
deed, various theories have postulated that AGN activity is 
responsible for truncating star formation (e.g., Springel et al. 
2005). While characterizing the star formation rates of broad- 
line AGNs is notoriously challenging, evidence for ongoing 
star formation has been seen in some cases (e.g., Canalizo & 
Stockton 2001). With narrow-line objects the galaxy is eas- 
ier to study directly, and Kauffmann et al. (2003) find a strong 
tendency for the most active objects in their sample of SDSS 
narrow-line AGNs to contain a significant post-starburst popu- 
lation. Ho (2005) use the [O II] doublet, which is intrinsically 
weak in high-ionization Seyfert galaxies, to find that star for- 
mation is suppressed in luminous AGNs, even when molecular 
gas is present (see also Kim et al. 2006); Ho speculates that 
this may be a concrete manifestation of AGN feedback. Our 
AGNs, while of low-mass and correspondingly low-luminosity, 
nevertheless have the highest growth rates of any local AGNs 
because of their high Eddington ratios (Greene & Ho 2007b). 
Unfortunately, we do not have a large amount of information 
about the stellar populations, both because of the AGN domi- 
nance and the low typical signal-to-noise ratio of the spectra. 

From the host galaxy colors, it is clear that recent star for- 
mation has occurred in some fraction of the galaxies. A more 
compelhng argument for intermediate-age stellar populations 
comes from the strength of the higher-order Balmer Unes ap- 
parent in the stacked spectra (Fig. 6). Balmer absorption is 
present in all Eddington ratio bins, but we cannot, from these 
data alone, quantify whether there is any correlation between 
AGN luminosity and strength of Balmer absorption. We are 
wary to measure stellar population indices for these objects 
without a clean measurement of the AGN luminosity, but the 
absorption-line spectra (particularly in panels c and d of Fig. 6) 
clearly contain a contribution from A stars. 

Each of these clues only places constraints on star forma- 
tion in the past ~ lO' yr. As always, it is far more challenging 
to find evidence for ongoing star formation, since the emission 
from the AGN itself generates nebular emission. The HST mor- 
phologies do suggest that star formation is ongoing at least in 
the spiral arms of the disk galaxies in the GH sample, but a 
more quantitative estimate is desirable. This was the motivation 
of Ho (2005) to use the [O II]/[0 III] ratio as an indicator of con- 
current star formation. Unfortunately, as we have seen above, 
we do not have a strong handle on the physical properties in 
the NLR, and therefore cannot uniquely interpret the observed 
spread in [O II]/[0 III] ratio. There is a mild correlation between 
the [O II]/[0 III] ratio and the ^ - r color (p = -0.3, P < 10"^; 
including or excluding the c sample), which one might expect if 
the [O n| arises from current star formation. On the other hand, 
if the NLR conditions correlate with Hubble type (not an unrea- 
sonable possibility), and hence indirectly with galaxy color, a 
similar correlation may result. With these data alone, we cannot 
place strong constraints on the ongoing star formation in these 
systems, nor are constraints on star formation alone sufficient 
to implicate a causal connection with bulge growth. We would 
like to determine whether there is an excess of star formation in 
the AGN hosts as compared to inactive dwarf galaxies. A com- 
bination of high spatial resolution optical imaging with HST 
and broader spectral modeling of the host colors should be a 
first step toward this goal. 
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Fig. 8. — Distributions of g-band absolute magnitudes. The magnitudes 
pertain to the AGN alone (top), the host galaxy alone (middle), and the en- 
tire system (bottom). Filled histograms represent objects below the detection 
threshold defined in Greene & Ho (2007b; the c sample). The AGN contribu- 
tion was estimated from Lhq. using the Lua-Lsim relation of Greene & Ho 
2005b, and assuming a spectral shape of /3 = —1.56 (fx oc A"''; Vanden Berk 
et al. 2001). Alternatively, the AGN magnitude may be estimated from the 
SDSS PSF magnitude, which leads to the distribution shown in dotted lines in 
the middle panel. 

4. SUMMARY AND FUTURE WORK 

We have increased the sample of low-mass BHs in galac- 
tic nuclei by an order of magnitude (174-1-55 additional candi- 
dates). In general the sample properties are completely consis- 
tent with that of the original GH sample. The objects are radi- 
ating at ~ one-third of their Eddington limit, and many of them 
have ROSAT counterparts. On the other hand, as expected, very 
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few have radio detections in the FIRST survey. 

Much like our previous sample, the current objects are found 
in sub-L* galaxies, suggesting that indeed low-mass BHs are 
found in low-mass stellar systems. The colors and composite 
spectral properties indicate that there are intermediate-age stel- 
lar populations, but we have no direct constraints on ongoing 
star formation. Although the SDSS imaging is neither deep 
enough nor of high enough angular resolution to determine the 
host galaxy structure, we do have relevant external data for the 
original GH sample. The bulge stellar velocity dispersions pre- 
sented in Earth et al. (2005) show that the objects are consistent 
with the low-mass extrapolation of the Mbh-o"* relation. There 
is a suggestion, from Greene & Ho (2006), that the slope of the 
AGN Mbh-ct* relation flattens somewhat at these low masses. 

While the systems (at least to first order) obey the Mbh - 
relation, they do not necessarily contain classical bulges. NGC 
4395 is an Sdm spiral galaxy (i.e., no bulge), while POX 52 is 
probably a dwarf spheroidal. As a result, the low-mass systems 
do not appear to obey the low-mass extrapolation of the Mbh- 
-Lbuige relation (e.g., Marconi & Hunt 2003). The total galaxy 
luminosities are certainly larger than one would expect, based 
on their BH masses, and this appears to be true even in the ab- 
sence of disk components (Earth et al. 2004; J. E. Greene et 
al. in preparation). Therefore even near their Eddington limits 
these BHs are not capable of providing as much energy per unit 
host galaxy mass as more massive systems. The mystery, then, 
is why they appear to obey the same Mbh - relation. 

Along these lines, we should note that the revised radius- 
luminosity relation slope (Bentz et al. 2006) substantially in- 
creases the virial masses at low BH mass, which in turn in- 
creases the apparent flattening in the Mbh - relation slope at 
low mass previously reported by Greene & Ho (2006). More 
insight should be gained into the reality of this slope change 
once we have measured for the present sample. Ultimately, 
we wish to elucidate the degree to which the slope and (cru- 
cially) intrinsic scatter of the Mbh - relation has any mass 
dependence (e.g., Robertson et al. 2006; Greene & Ho 2006). 



Indeed, there are a number of interesting applications to pur- 
sue with the improved statistics afforded by this new, larger 
sample. We plan to examine the environments and clustering 
properties of these objects, and see how they compare to in- 
active galaxies of similar mass and color. At the same time, 
we would like to constrain the stellar populations of the host 
galaxies in a more quantitative fashion. For the GH sample, 
at least, this will be enabled by a combination of HST imag- 
ing and Spitzer IRS spectra. Most immediately, we will use 
the HST imaging, in combination with stellar velocity disper- 
sion measurements, to place the GH sample on the fundamental 
plane. In combination, we will use these data to examine dif- 
ferences in the evolution of low-mass systems as compared to 
supermassive BHs in massive bulges. 
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APPENDIX 

VIRIAL BLACK HOLE MASS MEASUREMENTS 

The virial technique is quite indirect, and its calibration is a matter of ongoing research. Therefore, it is worth describing our 
methodology in some detail. In particular, we discuss the potential systematic errors incurred by uncertainties in the slope of the 
radius-luminosity relation on the one hand and the derivation of robust gas velocity dispersions on the other 

At the present time, the radius-luminosity relation is calibrated using the ^ 30 objects with BLR radii determined from reverber- 
ation mapping. Bentz et al. have used HST to determine the true AGN luminosity, as free as possible of host-galaxy contamination. 
Their quoted radius-luminosity relation slope is significantly shallower than that of Kaspi et al. (2005). It is worth keeping in mind 
the important caveat that the luminosities of the objects in that study barely overlap with those in this work, and most are significantly 
higher than those considered here. While programs to measure BLR radii for lower-luminosity sources are underway, at the moment 
we have no choice but to adopt the results for the more luminous sources. 

As described in detail in Greene & Ho (2005b), the line luminosity provides a more robust measure of the AGN luminosity in 
the presence of host galaxy contamination (or contamination from non-thermal jet emission). Therefore we have used the Lsioo-Lua 
relation from Greene & Ho (2005b) in combination with the revised radius-luminosity relation from Bentz et al. (and uncertainties 
from M. Bentz, private communication) to derive a final BH mass estimator: 

/ . X 0.45±0.03 /r7\T.TU\/i \ 2.06±0.06 

«»« = "»»x>»'(to^) (^) 

where we assume a spherical BLR (/ = 0.75). Here we are using the FWHM of Ha as a measure of the gas velocity dispersion. 
In principle the width of the broad emission lines represents the gas kinematics, but in practice it is unclear how to derive the most 
robust velocity dispersion measure. It has long been clear that the broad line shape depends systematically on other AGN properties 
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including line width (or BH mass) and luminosity (e.g., Sulentic et al. 2000); in general, narrower lines have more Lorentzian shapes 
(e.g., Veron-Cetty et al. 2001). If the geometry and kinematics of the BLR change with Eddington ratio (as suggested by these 
observations) then the assumption of a single geometrical factor / breaks down. In other words, the measured virial mass for a given 
BH may vary as changes in the Eddington ratio lead to changes in the BLR kinematics and thereby /. An extreme example of this 
effect is reported for AGNs with double-peaked Balmer lines (Bian et al. 2007). 

Without independent information it is impossible to determine / as a function of Lboi/^-Edd- Luckily, the Mbh-ct* relation provides 
an independent handle on the BH mass. Especially for the relatively low-luminosity AGNs considered at low redshift, it is safe to 
assume that the systems intrinsically obey the Mbh-o", relation, and adopt the Mbh-o", -derived mass as the true BH mass. Both 
Onken et al. (2004) and Greene & Ho (2006) find evidence for an average offsefl of -0.26 dex from masses derived with / = 0.75. 
However, as argued above, we expect that / is not a constant, but rather depends on Lboi/i^Edd- While Greene & Ho (2006) looked 
for evidence of a trend between / and Lboi/i'Edd, their dynamic range in both BH mass and luminosity was too limited to draw robust 
conclusions. 

Collin et al. (2006) perform an exhaustive study of the relation between line shape and virial mass. They use the Onken et al. 
sample of 14 AGNs with both reverberation-mapped masses and cr* measurements to investigate the merits of different measures of 
the line width. Specifically, they compare virial BH masses estimated from the FWHM and the luminosity-weighted second moment 
of the velocity profile (a) with the cr* -derived BH mass. As described above, they find that narrower lines are peaky (small values 
of FWHM/cr) while broader lines are boxy (large values of FWHM/cr). As a result, the inferred / may depend not only on whether 
FWHM or a is used to measure line width, but also on the line shape. Using the comparison between virial mass and cr*, they derive 
a different value of / for each subset of objects (9 peaky and 5 boxy) using both FWHM and a. They claim that while / is constant 
across line samples when a is used as the velocity dispersion indicator, the FWHM provides a biased measure of the virial mass 
because / is smaller for the boxy subset. However, investigation of their table reveals that the values of / derived from FWHM 
measurements are also consistent (within the quoted errors) with a constant value. Although the result may very well be real, it is 
certainly not statistically significant at this time. 

Moreover, we question the premise of Collin et al. that the best measure of line width is the one that minimizes changes in /. 
Rather, since the BLR geometry really does change with Lboi/i'Edd, one should calibrate different / values for different average line 
shapes, when sufficient statistics exist to do so in a robust fashion. Collin et al. present an additional piece of information that we 
believe favors use of the FWHM over cr. Because they are using reverberation-mapped AGNs, many epochs of observations are 
available. One may derive a so-called rms spectrum that isolates the variable part of the line. As argued in, e.g., Peterson et al. 
(2004), it is preferable to measure the width of the variable component of the line, which hopefully reflects the kinematics of gas at 
the measured radii. As it turns out, the FWHM measured from the mean and rms spectra are more self-consistent than the respective 
a measurements. There appears to be a very broad, non-variable component of the line which biases a more strongly than FWHM 
(perhaps an optically thin component; e.g.. Shields et al. 1995). While the very broad component cannot be removed in single-epoch 
observations, at least the FWHM is less biased by it. Finally, it is quite clear from Figure 1 of Collin et al. that the scatter in a 
measurements is considerably larger than in the FWHM. The same result was reported by Greene & Ho (2006). 

Based on all of these considerations, the following choices were made in calculating the BH masses presented herein: 

• We use the radius-luminosity relation reported by Bentz et al. (2006). Note that this differs from Greene & Ho (2007b). 

• We use the FWHM to determine the gas velocity dispersion. 

• Because we believe / is not a constant, but rather depends on Lboi/i^Edd, and in order to facilitate comparison with the 
literature, we do not apply an / correction but use the fiducial assumption of a spherical BLR (/ = 0.75; e.g., Netzer 1990). 

Obviously, virial masses remain highly uncertain, and potentially systematically biased. In addition to the many uncertainties high- 
lighted above, it is not yet clear the degree to which the BLR is a flattened structure, which means that inclination effects cannot 
be properly accounted for, even in a statistical sense. We are hopeful that an increased sample of reverberation-mapped sources 
(e.g., Kollatschny 2003), combined with an increase in the number of measured a^, measurements in AGNs, will decrease the out- 
standing systematic errors in these masses in the near future. 



Note that the Greene & Ho measurement, while including many more AGNs, is still dominated by the reverberation-mapped sources, which have considerably 
smaller systematic errors in their virial BH mass estimates. 
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